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Abstract 

We introduce a method to lower bound an entropy-based measure of genuine multipartite entanglement via nonlinear en- 
tanglement witnesses. We show that some of these bounds are tight and explicitly work out their connection to a framework 
of nonlinear witnesses that were published recently. Furthermore, we provide a detailed analysis of these lower bounds in the 
context of other possible bounds and measures. In exemplary cases, we show that only a few local measurements are necessary 
to determine these lower bounds. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.65.Fd, 03.65.Aa 
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I. INTRODUCTION 



Quantum entanglement is central to the field of quan- 
tum information theory. Due to its numerous applica- 
tions in upcoming quantum technology much research has 
been devoted to its understanding (for a recent overview 
consider Ref. 

Especially in systems comprised of many particles entan- 
glement provides numerous challenges and of course po- 
tential applications, such as building quantum computers 
(see Ref. |2|), performing quantum algorithms (the con- 
nection to multipartite entanglement is demonstrated in 
Ref. Q) and multi-party cryptography (see e.g. Ref. Q). 
Furthermore, the understanding of the behavior of com- 
plex systems seems to be closely linked to the understand- 
ing of multipartite entanglement manifestations, demon- 
strated by the connection to phase transitions and ion- 
ization in condensed matter systems (e.g. [Bj]), the prop- 
erties of ground states in relation to entanglement (as 
shown e.g. in Ref. [E 0), or potentially even biological 
systems (such as e.g. bird navigation @). 
In order to judge the relevance of entanglement in such 
systems it is crucial to not only detect its presence, but 
also quantify the amount. The structure of entangled 
states, especially in multipartite systems [9j], is very com- 
plex and the question whether a given state is entangled 
is even NP-hard [T3|. Thus, in general, it will not be 
possible to derive a computable measure of entanglement 
that reveals all entangled states to be entangled and dis- 
criminates between different entanglement classes. Fur- 
thermore, full information about the state of the system 
requires a number of measurements that grows exponen- 
tially in the size of the system. For the detection of en- 
tanglement in multipartite systems most researchers have 
therefore made it a primary goal to develop entanglement 



witnesses, which via a limited amount of local measure- 
ments can detect the presence of entanglement, even in 
complex systems (for an overview of multipartite entan- 
glement witnesses consider Ref. (TT| ) . 
The expectation value of witness-operators are usually 
expressed in terms of inequalities, which if violated show 
the presence of entanglement. Nonlinear witnesses (first 
introduced in Ref. [l2| see also early discussions in e.g. 
Ref. [lH) provide a generalization that is no longer a lin- 
ear function of density matrix elements, but a nonlinear 
one. Thus one cannot reformulate the criteria in terms 
of an expectation value of a hermitian operator (unless 
one considers coherent measurements on multiple copies 
of the state, which out of experimental infeasibility we do 
not discuss in our manuscript). We will henceforth refer 
to inequalities that involve nonlinear functions of density 
matrix elements as nonlinear entanglement witnesses. 
Recently some authors pointed out a connection between 
the possible amount of violation of these nonlinear in- 
equalities and quantification of entanglement in multi- 
partite systems (in Ref. [l4| and Ref. [IH). 
The aims of this paper are twofold. First to systemat- 
ically show the connection of numerous witnesses to a 
meaningful measure of genuine multipartite entanglement 
and second to use this established relation for the devel- 
opment of novel witnesses, which by construction give 
lower bounds on that measure. To that end we follow 
and generalize the approach from Ref. [l5[ . 
It turns out that only a small number of density matrix 
elements enters into our lower bounds, making the con- 
struction experimentally feasible even in larger systems 
of high dimensionality. 
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II. A MEASURE OF MULTIPARTITE ENTAN- 
GLEMENT AND ITS LOWER BOUNDS 

A. A measure of genuine multipartite entangle- 
ment (GME) 

The entropy of subsystems has often been used, in or- 
der to quantify entanglement contained in multipartite 
pure states (e.g. see P, Il6l - [l9j |). In this paper we will 
follow the definition first presented in Ref. [ly] and define 
a measure of GME for multipartite pure states as 



E m {\i>) := nun yJs L fa) = min ^2 (l - Tr(p 2 )) , 

(1) 

where Si (p 7 ) is the linear entropy of the reduced den- 
sity matrix of subsystem 7, i.e. p 7 := Tr 7 (|^)(?/;|). The 
minimum is taken over all possible reductions 7 (where 
the complement is denoted as 7) , which corresponds to a 
bipartite split into 7I7. 

As any proper measure of multipartite entanglement for 
pure states can be generalized to mixed states via a con- 
vex roof, i.e. 



E, 



t (p) := inf 
{Pi,\i>i)} 



5>£ m (l^)(^l) 



Due to its construction this measure fulfills almost all 
desirable properties one would expect from measures of 
GME (see Ref. [lU for details). Because computing all 
possible pure state decompositions of a density matrix is 
computationally impossible even if one is given the com- 
plete density matrix, we require lower bounds to be cal- 
culable for this expression. 

Also note that a lower bound on the linear entropy di- 
rectly leads to a lower bound on the Renyi 2-entropy 



(2) f21 

Sjj (pj) via the relation S R (p y ) 



io g2 (^¥^), 



which also provides one of the physical interpretations 
of this measure. The Renyi 2-entropy in itself is a lower 
bound to the von Neumann entropy <S(p 7 ) and the mutual 
information can be expressed as / 77 := S(p 1 ) + 5(p 7 ) — 
S(p) = 2S{p 1 ). Thus by our lower bound we gain a lower 
bound on the average minimal mutual information across 
all bipartitions of the pure states in the decomposition, 
minimized over all decompositions. 



B. Linear entropy and its convex roof 

The state vector of an n-partite qudit state can be ex- 
panded in terms of the computational basis 



d-i 



w = E 



|«1,«2, ' ' ' 



where a basis vector is denoted by ry = (11,^2, ■ ■ ■ ,i n ) S 
N® n .This vector notation will facilitate the upcoming 



derivations. A crucial element of the notation in this 
paper will be the permutation operator acting upon 
two vectors, exchanging vector components correspond- 
ing to the set of indices. E.g. the permutation operator 
■P{i,3}( 7 ?i) V2) will exchange the first and third component 
of the vector r\\ with the corresponding component of the 
vector 772, i.e. 

P {1 , 3} (01213, 30121) = (31113,00221). 

Using this notation one can write down a very simple 
expression for the linear entropy of a reduced state p 7 
(derivation see section 1 1 al in the appendix) 



(2) 



S L (Ay) = E Hi' 



where (^7) ^2) = (wiiVz)- For pure states we can of 
course find lower bounds on E m {\ijj) by lower bound- 
ing the linear entropy for all possible bipartitions. For 
mixed states we can then provide a lower bound for the 
convex roof E m (p). We now illustrate our method in one 
exemplary case and then continue to articulate the main 
theorem. 

Note that the linear entropy of subsystems has been 
widely used for lower bounding measures of entanglement 
due to the well known and simple structure of eq.©. 
None of the previous methods, however, work for lower 
bounding the inherently multipartite measure E m (p) , due 
to the additional minimization over all bipartitions in 
each decomposition element of the convex roof. 



C. W-states 

In order to demonstrate how our framework works let 
us start by deriving the explicit lower bound detect- 
ing the three-qubit W state \W) = ^(|001) + |010) + 
1 100)). For three-qubit states there are three bipartitions 
(1 1 23, 2 1 13, 3|12) and thus we have three linear entropies 
to look at in order to calculate E m (\i/j) (ifj\) , 



VSl(px) = 



(3) 



2 Vl c ooiCioo - C101C000I 2 + Icoiocioo - c lw c mo \ 2 + (•••), 
VSdp~2~)= (4) 
2\/|coioCioo - CnoCoool 2 + IcoioCooi - ConCoool 2 + (•■•): 
\fS~M= (5) 

E c »)l ? ')' 2A/IC001C100 - C101C000I 2 + IcoioCooi - c iic oo| 2 + (■■■)• 

77GN®" 



Now using yja 2 + b 2 > -j= (a + b) (which is a specific case 
of the incquality |A.5l in appendix l 1 bp and \a— b\ > |o| — \b\ 



2 



it is obvious that 



2(|c 00 icioo| - |cioic oo| + |c ioCioo| - |ciio£ooo|) E m (p) > V2(| <001]/o|100> I + |(001|p|010)| + [ <100|/>l010> [)- 

(6) ^-((010^1010) + (100|/9|100) + (001|p|001»- 

-\/2\/(101|/?|101)(000|p|000)- 

(7) V2V(110|/9|110)(000|p|000)- 
^^/(OlllplOllXOOOlplOOO) . (12) 



V2 



2(|cqiqCioo| - |ciiqC oo| + IcqiqCqqiI - [cqiiCoooJ) 

V2 



2(|c oiCioo| - Icioicoool + Icoiocooil - IcoiiCoool) 



V2 



(8) 



Then using \ab\ — i(a 2 + b 2 ) < we can add one negative 
term for each entropy and it will still be a lower bound, 
i.e. we add |c ioCooi I — -|(|coio| 2 + |c oi | 2 ) in the first lower 
bound, |ciooCooi| - ^(|cioo| 2 + |c oi | 2 ) in the second and 
Icoiociool - jdcoiol 2 + |cioo)| 2 in the third. Then we can 
use that mm{P-Ni,P-N 2 ,P-N 3 ] > P-N 1 -N 2 -N 3 
and end up with 

E m (\ijj){i>\) > V2(|c 00 iCioo| + |c oiCoio| + |ciooCoio|)- 
/2 

^-(|c io| 2 + |cioo| 2 + |c oi| 2 )- 

V^dcioicoool + Icnocoool + IcoiiCoool) ■ 

(9) 

Finally we can bound the convex roof using the following 
two relations 

inf Vk|c; iC ; 2 |> \{m\p\m)\, (io) 
inf Vpil44J < V(m\p\m)(m\p\m), (n) 

{pi,m)} *-r 

and end up with a lower bound for mixed states as 



Surprisingly this leads directly to the nonlinear entan- 
glement witness inequality presented in Refs. [13, Hl| up 
to a factor of Using only simple algebraic relations 
we have thus shown how to lower bound the convex roof 
construction. The first apparent strength of this lower 
bound is the limited number of density matrix elements 
needed to compute it. E.g. in our exemplary three-qubit 
case only ten out of possibly sixty-four elements need to 
be measured. Obviously we can extend the analysis using 
the same techniques to systems beyond three qubits. 



III. A GENERAL CONSTRUCTION OF LOWER 
BOUNDS ON THE GME MEASURE E m 

Now we can generalize the connection of the 3-qubit 
W state witness and the measure E m . Just as for three 
qubits we can always get lower bounds by summing the 
coefficient pairs c^c^ that belong to a certain target 
pure state and appear in some or all reduced linear en- 
tropics. The construction of such general lower bounds 
also starts by selecting a subset of coefficient pairs that 
will be translated into off-diagonal elements p Vl ,ri2i where 
(771 , V2) is the vector basis pair denoting the row and col- 
umn of the element in density matrix p. We denote the 
selected vector basis pairs as R := {(771,772)}- Then we 
can repeat the steps analogously to eq.(6-ll) and arrive 
at a general lower bound on the measure as the following 
theorem: 
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Theorem 1 (A general lower bound on the GME measure) For a set of row- column pairs R = {(?7i,?72)}, the 
genuine multipartite entanglement measure E m has the following lower bound: 



Em. > 2 




E ( \PviV2\ 



E yfav 

7Sr(r/i,T; 2 ) 



[\ E N n\Pw\ 



(13) 



The right- hand- side of eg. [TS\) defines a GME witness Wr(p), where p nitV2 := (r)i\p\ri2}, {rili'Hl) := ■^('Tij'Ta); 
r(?7i , 772) := {7 : (nj,!]^) (£ R} and I(R) '■= {n : 3rf that (rf ,rf) or (q,rj ! ) £ R} is the set of basis vectors n, which 
appear in the set R. 

Nr is the maximal (or minimal) value of |i? 7 | over all possible bipartitions 7I7, where R 1 is the set of coefficient 
pairs {c m ,c m ) € R, which do not contribute to the ^-subsystem entropy. 

N n are normalization constants given by the maximal value of over all possible bipartitions 7)7, where ri^ is the 
number of coefficients c n from some pairs in R, which are not counted in the ^-subsystem entropy (and how many are 
counted depends on whether one chooses Nr to be maximal or minimal). 



Proof. See Appendix [2] for the full proof. ■ 

It is evident that not every choice of coefficient pairs will 
yield a useful lower bound, because one really needs to 
select those that are actually contributing to multipar- 
tite entanglement. There is however always an obvious 
choice. The set of coefficient pairs R must be chosen such 
that in every subsystem at least one of the elements of 
R contribute to the linear entropy of the reduced state. 
E.g. in the case of GHZ states given in a specific ba- 
sis \GHZ) = ^(|0)® n + \\)® n ) one would choose the 
pair (00 • • • 0, 11 • • • 1), which contributes to all reduced 
entropies. In the general case however there is still some 
freedom of choice left to get a valid lower bound. For 
some sets R it can happen, that the coefficients do not 
contribute to every subsystem entropy equally (which we 
show in an exemplary case in section ft V A j) . Then one can 
choose Nr in different ways, but in all considered cases 
we found that choosing it maximal or minimal will pro- 
duce the best bounds (where choosing it maximal usually 
yields the tightest bounds close to pure states, whereas 
choosing it minimal improves the noise resistance) . Since 
these coefficients are in general basis dependent, so is also 

our witness construction. The prefactor \ R ^ N J sug- 
gests that the optimal basis for constructing such a lower 
bound is given by the minimal tensor rank representation 
of the pure state. 



IV. APPLICATIONS AND EXAMPLES 
A. Four-qubit singlet state 

Let us illustrate how to apply Theorem [T] with an ex- 
plicit example. In an experimental setting where one ex- 
pects to produce a four-qubit singlet state (which was 
e.g. discussed in the context of solving the liar detection 



I 

problem in Ref. [13]), i.e. 

\S 4 ) =-l=(2|0011) + 2|1100)- 10110) 

— |1001) — |1010) — |0101» , (14) 

one is confronted with the following expected coefficients: 

cooii,ciioo,coioi,cioio,c iio,ciooi- Following the recipe 
of theorem Q] we now select some coefficient pairs. Wc 
could choose e.g. i?i = (0011,0101), R 2 = (0011,1010), 
R 3 = (0011, 0110) and R 4 = (0011, 1001), such that R = 
{Ri, R2, R3, Ra}- For this selection we use theorem [T]to 
bound the GME measure. We see that in every subsystem 
at least two of these pairs appear naturally. Although 
there are more coefficient pairs we now choose to only 
take into account two per subsystem entropy and thus 
choose Nr to be the minimal number of coefficient pairs 
in every subsystem which gives Nr = 2. Thus we need 
to add negative terms that compensate for the missing 
terms just as we did in the three-qubit case, but now we 
need to do it two times in every subsystem. This results 
in the following individual prefactors N v for the diagonal 
elements: iVoon = 2 (as this coefficient appears in two 
missing pairs in every subsystem), iVoioi = 1, A1001 = E 
-^Vloio = 1 an d -A0110 = 1 (as those appear maximally 
once per subsystem entropy). Inserting this in theorem [T] 
we end up with the lower bound as 

2 

E m {p) > -y=(\p Rl \ + \PR 2 \ + \pR 3 \ + \PR 4 \ 

— v^Pon 1,011 1P0001, 0001 — V^ 0111 ' 0111 ^ 0010 . 0010 

— y/PlOll, 101 1^0001,0001 — \/Pl011. 101 1P0010, 0010 

— 2 (Poioi.oioi + Piooiaooi + Pioio,ioio + Pono.ono) 

— P0011.0011) • (15) 

We have thus created a nonlinear witness function that 
lower bounds our measure. From an experimental point 
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of view this is very favorable as few local measurement 
settings suffice to ascertain the needed thirteen density 
matrix elements (especially since the nine diagonal ele- 
ments can be constructed from a single measurement set- 
ting). Of course we could also exploit the connection of 
our lower bound to the Dicke state witness Q 2 (which 
is discussed in section HV C[) . which also detects GME in 
this state (although at the cost of more required mea- 
surements). In this case even the resistance to white 
noise is more favorable with our construction method, 
as for a state p = p\ S±) (S4I + ij^pl this exemplary lower 
bound detects GME until p = |^ ~ 0.72, whereas the 
old witness construction yields a worse resistance up to 
p = || w 0.77. This shows the versatility of our general 
approach. By choosing certain coefficients one can tailor 
these lower bounds to specific experimental situations. If 
one is confronted with a low noise system it is always 
beneficial to choose as few coefficients as possible, such 
that very few local measurements suffice (even a number 
that is linear in the size of the system is often sufficient). 
Every additional measurement can then be included in 
the lower bound and improves the bound and its noise 
resistance if necessary. 

B. Bipartite witnesses and lower bounds on the 
measure 

Although we have presented our theorem and measures 
in the general case of n-qudits, we can always apply the 
lower bounds also for n = 2, as our theorem holds for 
any n and d. Suppose we are given a bipartite qutrit 
system and want to lower bound the concurrence with 
only a few local measurements. If the expected state is 
e.g. \ip) = -i=(|00) + |11) + |22)) wc can use the lower 
bounding procedure outlined above, yielding 

E m {p) > -|(3?e[(00|p|ll)] - V(01|p|01)(|10 P |10) + 

Ke[(00|/o|22)] - V <02|p|02) (20|p|20) + 

&e[(ll|p|22)] - V<12|p|12) <21|p|21>) . 

(16) 

In order to determine the lower bound we have to mea- 
sure nine different density matrix elements. Of course 
any density matrix element can always be obtained via 
local measurements. How these measurements can be 
performed in a basis consisting of a tensor product of the 
generalized Gcll-Mann matrices we show explicitly in ap- 
pendix [31 

It turns out that these nine different density matrix el- 
ements can be obtained via ten local measurement set- 
tings. Let us study the lower bound in the presence of 
noise. Suppose we have white noise in the system, i.e. 
P = p\4'){4 ! \ + "ipl- Calculating the lower bound results 
in E m (p) > 2 ^ 4 ^1 1 ^ , which is equivalent to the analytical 



expression of Wootter's concurrence for these systems (as 
proven in Rcf . [H, HH ) . In this case we have a necessary 
and sufficient entanglement criterion and a tight lower 
bound on the concurrence from ten local measurements 
for a special class of states. Indeed if one generalizes this 
example to arbitrary dimension d, we find that the bound 
is always tight for bipartite isotropic states. 



C. Dicke States 

We will now continue to show how this construc- 
tion relates to an entanglement witness for Dicke-state, 
which arc multi-dimensional generalizations of the W 
states (which were first introduced in the context of laser 
emission in Ref. [H[). 

In the original article fl5l |. where this approach was 
first introduced, the authors connected the violation of 
a witness suitable for GHZ states (first introduced in 
Ref. [13] and later presented in a more general frame- 
work in Ref. [Hj]) with a lower bound on the measure 
E m . We want to follow this approach and establish a 
general connection between a set of witnesses suitable for 
all generalized Dicke states introduced in Ref. [H[ and 
generalized in Ref. [H| . To that end let us first introduce 
a concise notation for those states. 

Let a be a set containing specific subsystems of a multi- 
partite state. We then define the state \a l ) as a tensor 
product of states |Z) for all subsystems not contained in a 
and excited states \l + 1) in the subsystems contained in 
a. E.g. for the four-partite state |{1, 3} 2 ) we have 1 3232) . 
Using this abbreviated notation wc can define a gener- 
alized set of Dicke states, consisting of n d-dimensional 
subsystems, as 

\D d j= 1 E E i ai >< 

V(m)( d_1 ) l=0<x:\ a \=m 

where the parameter m denotes the number of excita- 
tions, with < m < n. 

Since the explicit form of the nonlinear witness from 
Ref. [H| will be used in the following considerations we 
will repeat it in appendix [H For all biseparable states 
this witness Qm is strictly smaller equal zero, i.e. 

Q (p) < <= p is biseparable 

Q (p) > =>■ p is multipartite entangled . 

Furthermore, the witness can also detect the "dimen- 
sionality" of GME, by which we mean the maximal num- 
ber of degrees of freedom f p (f p < d) that occurs in the 
pure states of an ensemble constituting p, minimised over 
all ensembles (this is the natural generalization of the 
concept of Schmidt number 12511 to multipartite systems, 
further explored e.g. in Rcf. )24| ) . I.e. the dimensionality 
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D. PPT- Witness and Our Witness 



FIG. 1: (Color online) The witness Q^? can quantify the gen- 
uine multipartite entanglement with the GME dimensional- 
ity f p , where f p < d . A state is biseparable iff f p = 1 
and GM entangled iff f p > 2. For Dicke states it holds 

Q%>(\DL){DL\)=f-h 



inf max(min(rank(/5 7 ))) 



(18) 



is defined as 
Since 

Ql? (p)</ P -l,V P , 

we can directly infer that 

(p) > / - 2 => / p > / 

In fig|T] we show how Qm detects the GME dimen- 
sionality. The maximal violation of these inequalities 
is always achieved for m-excitation Dicke states, i.e. 
Q%\\Di)(Di\) = d-l. 

If we can find a proper R, as a result of theorem [T] 
that uses the Dicke state coefficients, we can connect a 
lower bound of the measure E m with the GME witness 
Qrn (p)- Indeed choosing the ordered subset R a of the 
set of coefficients a used in (|A.23[) , i.e. 

R a = {(a Q ,/3 b ) e a : a < b} , 

we immediately arrive at a lower bound on E m as 



E m (p) > rriA 



\Rtr \ — Nk 



Qln (P) > m s 



(19) 

where |i? CT | = ^ (d — 1) (™Jm(n — m). In this case iV I( < 
m(n — m — l)+0 (d — 3) (n — m), where O is a Heavisidc 
step function. 




W R (p) - undetectable states 



(p\Xi)(Zi\ T ') = 



FIG. 2: (Color online) PPT witness compared to our lower 
bounds (given in terms of the nonlinear witness Wr (p)): The 
set of Wr (p) undetectable states denotes the set that is not 
detected by one specific Wr(p) and is strictly larger than 
the set of PPT-states. However the set of states detected by 
Wr (p) is strictly larger than the set detected by any standard 
PPT-witness. 



Using the result on entanglement across bipartitions 
from the previous section we can explore the relation of 
our lower bounds to other bipartite entanglement wit- 
nesses. In our witness construction, the permutation op- 
erator P-y acting on a pure state is a 7 ^-partial trans- 
pose operator, i.e. P 7 |i/')(^/ , | = (\i>) (ip\) T ~*^ (in the sense 
that our permutation operator now acts upon the index 
pairs of the coefficients of the pure state). It is thus 
intuitive to believe that there is certain connection be- 
tween our witness and a PPT-witness [26j]. Indeed our 
witnesses are related to a standard PPT-witness con- 
struction (where the witnesses separate the convex set of 
states that are positive under partial transpose (PPT) 
from its complement). E.g. for diagonal GHZ states 
we can use the standard PPT-witness construction which 
goes as follows. For |GHZ^ l!??2 ) := ^= (\r)i) + (772)) with 

Vi + ? /2 = {d — 1, ■ ■ ■ , d — 1), we can use the eigenvector 
belonging to the negative eigenvalue of the 7|7-partial 



transposed |GrIZ, n ^ 2 ) (GHZ m ,, )2 1 7 which we denote as 



A: 



1 

'm,V2/ ~ v/2 

PPT-witness and write its expectation value as 



,> 



(| — \V2 ))■ One can then construct the 



«pS(p>I^ 1 ,, 2 »=^(|a-,, 2 >(a 



(20) 



For instance in the thrce-qubit case, 



\\~ \ /\~ l Tl ' 23 - 

I A 001,110/ \ A 001,110 I 



G 



/ 




\ 

-1 



Vo 



-1 



(21) 



0/ 



With the PPT- witness construction in ea. (|2TJ]) we end up 
with the following PPT-witness expectation value 



n 



~t\i 
ppt 



Pi \ X GHz)) 



1 



'Vi Vi 



P 



Vi Vi 



Re(p mm ) 



(22) 



Under the fixed bipartition 7I7, we construct our witness 
by choosing R = (171,772) as 



W? , 

(VI, V2) 



(p) = y/PnlvJP 



\Pvi' 



(23) 



It is obvious that 



-Kuv.) 



(p) 



< 



^ppt (Pi l^GHZ 



Hence we say that the witness Wr (p) is stronger than 
the PPT-witness fig (p, |A" )I72 ». 

The relation between our witness, the PPT-witness and 
the PPT-convex set is illustrated in figlH For clearness 
we just draw two PPT-witnesses in the figure. For the 
n-qudit case there are \d n such eigenvectors |Ar ), 
corresponding to negative eigenvalues. Every witness 



O 



7I7 
ppt 



p, \\ is tangent to the set of PPT states (i.e. 



there exists one PPT state for which the witness yields 
zero). However also our witness Wr (p) is zero for all 



these PPT states, i.e. our new witness detects more states 
than the traditional PPT-witness. 



V. CONCLUSIONS 

In conclusion we have presented a method to derive 
lower bounds on a measure of genuine multipartite en- 
tanglement. We show that in experimentally plausible 
scenarios (i.e. one knows which state one aims to pro- 
duce) we can derive such lower bounds simply based on 
coefficients of the corresponding pure states. We also con- 
nected the lower bound construction to a framework of 
nonlinear entanglement witnesses developed in Refs. [20I — 
|24| . These witnesses are experimentally feasible in terms 
of required local measurement settings. We provide fur- 
ther evidence in the bipartite case, where we also show 
that for certain families of mixed states our lower bounds 
are tight. 

Some open questions remain, such as whether this gen- 
eral construction method will work for all kinds of states 
and how it can be generalized beyond just multi- and bi- 
partite entanglement, but anything in between. We want 
to point out that recently also other authors have used a 
similar approach to bound this measure in the bipartite 
case and for multipartite W states j30| . 
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Appendix: Proofs 



b. An Important Inequality 



1. The Formulas Used in The Main Article 

a. Reduced linear entropy of pure states 

Let = X^eN®" °n \ v) De an n-qudit pure state. The 
linear entropy of \rp) can be written as 



S l (p 7 ) = £ 



C Vl Cr) 2 C^t- C^t- 



(A.l) 



where (rf[,T}Z) = P 7 (771,572). 

Proof. The linear entropy regarding a specific parti- 
tion 7I7 is defined as Sl (p 7 ) = 2(1 — ir (p 7 )), where 
p 7 is the 7-reduced matrix of p. The trace of p 7 is 
* r (pV) = „ (p-y) (p-v) , where i/^ is the 
subspace of the reduction 7. We separate the summation 
into diagonal and off-diagonal parts. For the diagonal 
part we use the normalization condition to evaluate its 
value. 



£ (P7) Ql a 2 (P7) a 2 ai 
: 1- £ l C "l<»/3l| 2 |Ca 2 «./3 2 | 2 

c ctl <S>/?1 c a 2 ®/3l Ca 2 ®/3 2 C ai (g>,8 2 ■ 

ai^a 2 S-f/-y 
/3l,/3 2 £ff, 



(A.2) 



By exchanging the indices ol\ and 012 one has 

t r (P 7 ) = ^ — 9 £ l c ai«>/3i c a 2 «i^2 — c ai®/3 2 c ci2®/9il 



= 1 - 



ai#a 2 £-H 7 
/3i,/3 2 eff, 

- r 



(A.3) 



where r] = a <g> ft and (777,772) = -F7 (771,772). The linear 
entropy is then calculated to 



S L (P 7 ) = 51 

1715^172, en d n 



c m c rj 2 c r)1 c r)1 



(A.4) 



The following is an inequality, which is crucial for 
derivation of the prefactor W r^rrjVfl nl ^ ne theorem [TJ 



i/i£n 2 > 



(A.5) 



Proof. We prove this inequality by constructing two 
vectors as follows (using |/| = n) 



( «i \ 




> n times 



,y 



J 



(A.6) 



V < J 



The right hand side of IA.5I can be written as the scalar 
product of x and y. 



= Y,a ia * = \x-y\. (A.7) 



According to the Cauchy-Schwarz inequality, one can de- 
rive 



(A.8) 



2. Proof of Theorem [T] and Approach of Construc- 
tion of a GME Witness 

Firstly one can estimate the lower bound on Sl (p 7 ) 
by summing its elements over a selected Region R, and 
dropping the other non-negative summands (i.e. lower 
bounding them with 0), 



Sl(p 7 )>4 £ 14,4 -c; 

-4 £ 

(m^2)e-R\-RT 



(A.9) 



Here we add a prefactor 4 in eq. (|A.9[) , since the symmet- 
ric factor of all (771,772) equals 4. That means for every 
(771,772) there are three other (771,772) having the same 
value of \cf ll Cfi 2 — c^c^l as (771,772). Here we choose a 
non-degenerate vector basis set i?, and therefore need a 
prefactor 4 in the lower bound. The set W is the subset of 
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R, whose elements do not contribute to the linear entropy, 
i.e. i? 7 := {(771,773) G R : (vj,^) = (VuVs) or (r) 2 ,r]i)}. 
Now we use the inequality (IA.5I) to bound the square root 
ofSi(p*). 



St 




P c 1 c L \ 



(A.10) 



Ry\ ^ I 



(A.ll) 



According to eo. pi A|) together with ea. (|A.ll[) . the lower 
bound reads 



2 inf y Pi 

{w.lV'i)}^ 



1 



r\ - \Ry* 



■ E 

V1.V2SR 









(\ C V1 C V2 1 




) 









(A.12) 

where 7, is the partition in which the linear entropy 
(jV'i) (V'il-y) °f IV'i) (V'il nas its minimum. By defin- 
ing the normalization factor Nr := min 7 |i? 7 | , which is 
the minimal value of |i? 7 | over all possible bipartitions 
{7I7}, we can extract the prefactor from the convex roof 
summation. 



1 



inf 

\R\-N R { Pi ,Wi> 



E 

VI >V2^R 



(A.13) 



The most difficult part of detecting entanglement of 
mixed states is a result of the mixing of the decomposi- 
tion coefficients c^c^. In the lab we have only the in- 
formation about the mixed density matrix element p Vim 
but not cL c l V2 , therefore we must exchange the two sum- 
mations in eq. (|A.13p . and mix the coefficients c^c^ into 
density matrix elements. Therefore we estimate the sum- 
mands with a bound, which is independent of the specific 
partition 7i|7;, by adding a summation of non-positive 

into the sum- 



terms Y^R-n 
mands. 

E 

(vi,V2)eR\R">i 



> E 

(vi,V2)£R 



'71 I 



n- 2 



"m 1 



4E(I4 

R~ 

> E 



i |2 



- E I 

7er(?7i ,172) 

,2> 



E 



I 1)1 "72 I ^ I ">7i _ ») 2 ' 

(V1,V2)GR \ 76r(77i,i7 2 ) 



i E <K! 2 : 



vei(R) 



(A. 14) 



where /(i?) := {?/ G Nf" : 3 (77, t/') or (77', 77) G R} is 
the set of indices contained in the set R, r (771 , 772) = 
{7I J 3 (771 , 772) (I R} and n 7i is the number of vector pairs 
in R li containing index ?y. In order to eliminate the de- 
pendence of the partition 7* , we define the maximal value 
of ?i 7 over all possible partitions {7I7} as N v := max, 71 7 . 



Then one can estimate the GME measure with eq. (|A.13l 
and IA.14[) as 



E m {p) 



>2 



1 



\R\ - N R 



E 

VI >V2^R 



inf >^ Pi 



E 

7er(i7!,?72) 



1 E 

v£I(R) 



(A.15) 



Now one can safely exchange the summation in eq. (|A.15I) 
and lower bound it with the triangle inequality (i.e. 
Epi -Pi I c ?n c r? 2 1 — IPm'jal) an d the Cauchy-Schwarz in- 

< y/PvlvlPnlnl)- Fin ally 



equality (i.e. J2 Pi P 
we arrive at the result 



Si V2 



E m (p) 



>2 



R\-N F 



E 

V1.V2&R- 



PV1V2 



E • 

7er(i7i,77 2 ) 



E N vPvv 

riel(R) 



(A.16) 



where p mm := <?7i |p|-??2> ■ 

Above is the proof of theorem [T] in the case of Nr '■= 
min 7 |i? 7 |. For the choice of Nr := max 7 |i? 7 |, one 
just needs to calculate max 7 |i? 7 | at the first step, i.e. 
cq. (|A.9p . then pick up \R\ — max 7 |i? 7 | elements from 
i?\i? 7 as summation region in the second line and then 
repeat the whole proof above. At the end we will at- 
tain the same expression for the lower bound on E m as 
eq. (|A.16[) . but with different N n from the ones before 
Nr = min 7 |i? 7 |. N^ in this maximum choice is greater 
or equal to the one derived in the minimal-case. In the 
four-qubit singlet example in sec II V Al the value of N n is 
exactly the same for both choices. Therefore we choose 
the maximum, i.e. Nr = 2, to get a tighter lower bound 
on E m . 



3. Explicit decomposition of the bipartite witness 
into local observables 

The measurements needed to ascertain the relevant 
density matrix elements in the bipartite scenario can be 
performed in a basis consisting of a tensor product of 
the generalized Gell-Mann matrices. We continue to pro- 
vide for each of the density matrix elements above their 
respective coefficients. The density matrix elements are 
either off diagonal elements or diagonal elements. The 
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off-diagonal elements can be obtained by expectation val- 
ues of the symmetric and antisymmetric generalized Gell- 
Mann matrices: 



For further reference the coefficients are given as: 



A 



A 12 = 



1 o\ 

1 o o ,A 

0, 



1 

10 



/0-i0\ . 13 /0 0-i\ 

[ i o o ,A„ =oo o 
Vooo/ Vioo/ 



A 



A 23 = 



They can be written as follows: 



Me [<00|p|ll>] = - (Al 2 ® A s 12 - A* 2 ® A* 2 ) 



Ke[(00|p|22)] = -<A 
Re[<ll|p|22)] = i(A 



-Af-A^A* 3 ) 
A 23 -A 23 ® A 23 ) 



/ 
1 
V0 1 

(A.17) 


-i 
-t 

(A.18) 



(A.19) 
(A.20) 
(A.21) 



We now consider the terms obtained via the diag 
onal generalized Gell-Mann matrices. A° = ( o l o 



Ai = 



10 
0-10 



A 2 _ 



1 
10 
0-2 



1, 

We will expand the 



l 

,ooo/' " VS 
soughtafter terms into coefficients, utilizing the following 

basis: 

A2®A 2 

a|®a° 

A|®AJ 

Va^a 2 / 



6 = 



(A.22) 



(01 
(10 
(02 
(20 
(12 
(21 



|01) 
|10) 
|02) 
[20) 
|12) 
[21) 



b * 
b * 
b * 
b * 
b * 
b * 



1 1 



1 1 



1 



1 



1 



6'6V3' 6' 4' 8V3'6V3' 
11111 1 

11 11 



11 1 

l l 



3^3' 2^3' 
1 



r,0, 



6' 6^ 
3 



,0,0,0,- 



4^' 6^ 
1 1 



»0,-r 



4,o, 1 



3\/3' 2^3' 
1 



3V3' 6' '4^3' 6^ 



,0,-r 



4^3 
1 

4^ 

1 

6 
1 
6 
1 
6 



1 

'12 
1 

'12 



id) 



4. Explicit form of the GME witness Q 



Here we recall the explicit form of the nonlinear wit- 
ness from Ref. [23] . Using the notation for Dicke states 
introduced in section llV Cl we arrive at the following lower 
bound 



Q 



with 



m 



d-2 



EE <*W 



Zi'=0 CT 



SeA ) 1=0 a 



(A.23) 



m e {1, • • ■ , |n/2j} , iVo = (d - 1) m (n - m - 1) . 
a := {(a, (3) : \af\ (3\ = m - 1} , 

a ,V = l 

A:= I ,V <l 

\8\5c0\a) ,l'>l 



(A.24) 



r 



The properties of this witness are discussed in the main 
text. 
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